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Abstract 

We describe a novel use of the Dalitz plot to probe CP symmetry in three-body modes of 
B and D mesons. It is based on an observable inspired by astronomers' practice, namely 
the significance in the difference between corresponding Dalitz plot bins. It provides 
a model independent mapping of local CP asymmetries. We illustrate the method for 
probing CP symmetry in the two complementary cases of B and D decays: in the former 
sizable or even large effects can be expected, yet have to be differentiated against leading 
Standard Model contributions, while in the latter one cannot count on sizable effects, yet 
has to deal with much less Standard Model background. 
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1 Prologue 

While the announcement of the 2008 Nobel Prize in Physics has made official the status 
of KM theory as the main source of CP violation as observed in K and B decays, it does 
not close the chapter on it for three main reasons [1]: 

• We know baryogenesis in our Universe requires New Physics with CP violation. 

• A host of largely theoretical arguments suggests - persuasively in our view - that 
New Physics exists around the TeV scale with rich dynamical structures. In general 
those can provide several novel sources of CP violation. In that context one uses the 
high sensitivity of CP studies as a tool to search for New Physics and hopefully infer 
some of its salient features. One should keep in mind that observable CP asym- 
metries can be linear in a New Physics amplitude with the Standard Model (SM) 
providing the other one; therefore one achieves sensitivity to small contributions. 

• The LHCb experiment [2] is poised to acquire large sets of high quality data on the 
decays of B and D mesons. 

The stage for CP studies has recently become wider with the observation of D'^ — 
oscillations [3, 4]. Fortunately we can expect to continue our quest for CP violation with 
the continuing work of the Belle Collaboration, the hoped for realization of a Super-B 
Factory [5] and in particular with the beginning of the LHCb experiment. 

We will focus on three-body final states in the decays of B and D mesons and a novel 
strategy to probe CP symmetry in their Dalitz plots. Those two classes of transitions 
offer complementary challenges both on the experimental and theoretical side. 

• KM dynamics is expected to generate large CP asymmetries - to the tune of, say, 
10 - 20 % - in modes like B —>■ tttttt, Kttti, Kpp. Since those command very small 
branching ratios, wc still need very large samples of B mesons, as will be produced 
by the LHCb experiment. While we can be confident that such effects will be 
observed (supporting the experimentalists' enthusiasm), the real challenge at that 
time will be how to interpret a signal: does it reveal the presence of New Physics 
or is it consistent with being from KM dynamics alone? This will represent a non- 
trivial conundrum, since all the available evidence points to New Physics mustering 
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no more than non-leading contributions in B decays; this is often referred to as the 
" Flavour Problem of New Physics" . We will have a realistic chance to answer this 
question only if we have accurate as well as comprehensive data. 

On the other hand in most charm decays experimental bounds tell us we can hope 
for at best moderate size asymmetries even in the presence of New Physics. The 
redeeming feature is that the SM can generate CP violation at most at the 10~^ 
level in the 'best' cases and even significantly less in others. However we view the 
latter as good news: for almost any asymmetry observed in the near future would 
represent strong evidence for the intervention of New Physics. Using the language of 
'signal-to- noise' familiar from experimental studies: while we expect significantly 
smaller asymmetries in D than in B decays, the theoretical 'noise' or 'background' 
- i.e. contributions from the SM - is even more reduced in the former than the 
latter. Thus we conjecture 
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CP asymmetry 
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JJ decays 



B decays 

The charm branching ratios are typically sizable and we have already acquired a 
great deal of experience in describing them and their Dahtz plots. The central 
challenge then is to control systematics to a degree that allows probing asymmetries 
down to 10~^ or even better. 

So far no CP asymmetry has been established on the five sigma significance level. 
However we expect that to change soon and actually predict Dalitz studies to become one 
of the central tools for CP probes. 

• Accurate measurements of CP asymmetries will be a necessary (though probably not 
sufficient) condition for deciding whether they reveal the intervention of New Physics 
or not. In the case of charm decays the anticipated small size of effects constitutes 
the main challenge; in B decays, on the other hand, we have to learn how to subtract 
the presumably leading SM contribution. Dalitz plot descriptions with their many 
correlations yield o?;erconstraints, providing reliable validation tools. Tracking the 
time evolution in D^, Bd and Bs transitions characteristic of oscillations will further 
illuminate the underlying dynamics. 

• Establishing the intervention of New Physics in CP studies will of course represent 
a seminal achievement, yet we want to do even better. For our goal has to be to 
infer salient features of that anticipated New Physics. Some of those can be read 
off the flavour structure of the final states. Yet to infer the Lorentz structure of the 
underlying operator we have to go beyond final states consisting of two pseudoscalar 
or one pseudoscalar and one vector meson; for their amplitudes are described by a 
single number. However, once we have three pseudoscalar mesons or a baryon- 
antibaryon pair plus a meson, the kinematics is no longer trivial, and final state 
distributions can tell us whether spin-zero or spin-one couplings are involved in the 
transition operator. This important feature will be illustrated below. 
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Fortunately pioneering work has been done by Belle and BaBar: Based on a Dalitz plot 
study they have extended their probe of CP invariance to quasi-two-body channels in- 
volving resonances and indeed found intriguing evidence for a direct CP asymmetry in 
the mode i?^ — > [6, 7, 8], which - if established - would be a first. 

Relying on mass projections is only one way to use the dynamical information con- 
tained in a Dalitz plot, and it cannot be expected to harness the full potential of Dalitz 
studies. It is the interference between two neighbouring resonances that presumably pro- 
vides the most sensitive CP probe. For a direct CP asymmetry to surface one needs the 
interplay of a weak and a strong phase with the former in contrast to the latter changing 
signs under a CP transformation. For the latter one usually takes the strong phase shifts, 
which cannot be calculated from first principles. Yet when one deals with a finite-width 
resonance, its Breit-Wigner parameters can provide the required strong phase, which 
varies with the mass bin in a characteristic and largely predictable way. This feature can 
provide further validation for the experimental findings. Alas - a full-fledged Dalitz plot 
description requires huge statistics and considerable theoretical 'overhead' in selecting the 
transitions deemed relevant and parametrizing their amplitudes. It has to be the ulti- 
mate goal to develop such a complete description with as much accuracy as possible, yet 
that will be a long term task, and it is not clear what irreducible model dependence will 
remain. A full Dalitz plot description would help in the extraction of the CKM phases 
[9, 10]. There are, however, some theoretical issues concerning the large phase space of 
the B decay into three hght mesons [11] that need to be understood. 

To avoid such model dependence one can divide the Dalitz plot into bins, and then 
directly compare the CP conjugate Dalitz plot regions in a bin-by-bin basis. Yet results 
based on studying the ratio between the difference over the sum of the populations are 
quite vulnerable to fake effects from statistical fluctuations. Therefore we suggest a re- 
finement of such a direct comparison, namely to study the significance of the difference. 
This proposal has been inspired by what has become standard routine in astronomy when 
analyzing light sources in the sky [12]. Its main values lies in three aspects: 

• As illustrated later it provides a model-independent and robust method to determine 
CP asymmetries already with limited statistics and identify the regions of a Dalitz 
plot, where they occur. 

• This is particularly important when dealing with small or even tiny effects as ex- 
pected in charm decays. 

• Its findings provide powerful constraints on any full Dalitz plot model to emerge. 

In talking about 'limited' statistics we do not mean small statistics - a situation addressed 
in [13]. Since our method involves analyzing distributions even in sub-domains of the 
Dalitz plot, it requires substantial data sets. 

It has been estimated that LHCb will collect very sizable data sets of three-body 
decays already in one nominal LHC year: 

• about 10^ singly Cabibbo suppressed tt^tt+tt" and 10^ doubly Cabibbo 
suppressed K^n+n- /K^K+R-; 
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• around 10^ S± ^ K^n+n', K^K+K' and S± ^ tt^tt+tt-; 

• more than 10^ B"^ — > 7r^ir"'"X~ and the baryonic modes B"^ — > K^pp, Tr'^pp. 

The paper will be organized as follows: in Sect. 2 we will briefly review the basics of 
Dalitz plot analyses, introduce a novel observable for probing CP symmetry there and 
comment on the isobar model; in Sect. 3 we apply it to B mesons decaying into three 
light mesons and present Monte Carlo studies; in Sect. 4 we discuss analogous D decays 
where one has to face rather different challenges; in Sect. 5 we present expectations about 
in which direction and to which degree relevant theoretical tools might get refined in the 
near future before summarizing and giving an outlook in Sect. 6. 

2 Basics and Virtues of the Dalitz Plot 

2.1 Basics 

It is of course a mathematical triviality that local asymmetries are bound to be larger than 
fully integrated ones. Yet a Dalitz plot description translates such a general qualitative 
statement into a much more concrete one. For it exhibits all that can be learnt directly 
from the data on final states of three stable particles and their dynamics. Since the phase 
space density of the Dalitz plot is constant, any observed structure reflects the dynam- 
ics of the decay. Enhanced populations in certain mass regions can reveal the presence 
of a strong resonance and indicate their widths. The angular distributions characteris- 
tic for the spin of the resonances modulate the mass bands. Distorted or twisted mass 
bands point to the interference between resonances. These observations can be cast into 
a quantitative treatment by making an ansatz for the flnal state amplitude consisting of 
terms describing the moduli and complex phases of the contributing resonances and the 
non-resonant contribution. These entities contain a great deal of subtle dynamical infor- 
mation. Comparing them for CP conjugate transitions provides a very powerful probe of 
CP invariance. While CP violation has to enter through complex phases on the funda- 
mental level of the underlying dynamics, it can manifest itself in the Dalitz plot through 
differences in both the aforementioned moduli of the hadronic resonances and their phases 
for conjugate transitions. Since there are typically several resonances contributing to a 
decay, ample opportunities arise for CP violation to surface in a Dalitz plot. Hadronic 
'complexities' thus represent good news for the observability of CP asymmetries. They 
become a challenge only, when one undertakes to interpret a signal in a quantitative way. 
Yet even there a Dalitz analysis provides essential assistance: the reliability of Dalitz 
plot parametrizations can be inferred from the amount of overeonstraints they manage to 
satisfy. 

However Dalitz studies still retain a measure of model dependance due to the choices 
one makes concerning the resonances to be included and their parameterization and also 
due to the treatment of the non-resonant contribution; the S-wavc is the largest source 
of systematics due to strong dynamics. The greatly different phase space available in B 
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and D decays makes for an almost qualitative difference in how to treat them. We will 
comment on it later. 

While we maintain that such model dependencies can be reduced considerably with 
increasing data sets and, more important, with future theoretical insights, we want to 
propose a novel method for searching for CP asymmetries in three-body final states that 
is robust in two respects: it requires no model assumptions and provides an effective filter 
against effects due to statistical fluctuations. Yet flrst we will make a few rather technical 
remarks on how various phases enter the interference between neighbouring resonances. 



2.2 Phases with Breit-Wigner Resonances 

Due to CPT invariance CP violation can express itself only via a complex and pre- 
sumably weak phase. For it to become observable, we need the interference between 
two different, yet still coherent amplitudes. Oscillations can provide such a scenario - 
as can hadronization in general. The latter case is usually expressed by stating that 
the two amplitudes have to exhibit different weak as well as strong phases: Ai = 
ei^Te^^t^f)\Mi\ + e''*'Te'^2^f^\M2l with (0r>^f(/)) ^ (0r>^2*(/))- We can write it 
also in terms of phases that combine the weak and strong phases: 

M = e'^i(^)|Mi|+e'*^(^)|M2| , = 5f (/) + (2) 

M = e'^"i(^)|Mi| + (^~^^^f^\M2\ , Uf) = Sfif) - (t)^ (3) 

It is often implied that the strong phases Sf{f) carry a fixed value for a given final state 
/. This does not need to be true. More specifically it will definitely not hold when the 
final state contains a resonance. The Breit-Wigner excitation curve for a resonance it! 
reads 

^l^(s) = —9 TT^ > (4) 

introducing a sizable phase as expressed through 

ImFl^ (.) = , , "T/^^'^ , , (5) 
^ ^ ^ {ml- sf + {mRTR{s)Y ' ^ ' 

where r/j(s) denotes the energy dependent relativistic width. In our discussion of 
K^-K^-K^ we will focus on the interference between p° and /o to generate a CP asymmetry. 
The relevant amplitude components for and B~ are: 

M- = a^.e^'-Ff"^ cose ^aW-Ff^ (7) 

The 5± contain both the fixed weak and the strong phases with the Breit-Wigner functions 
^Bw introducing additional mass dependent strong phases as sketched above. For the 
/o we have followed BaBar's treatment [6] using the Flatte representation, which refiects 
the proximity of the KK threshold and the ensuing distortion of the resonance curve. In 
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Eqs.(6,7) above 6 is the angle between the 7r~ and the momenta, measured in the p 
rest frame. This angle describes the angular distribution of a vector meson. After taking 
the modulus square of these amplitudes one reads off that a CP asymmetry will arise, 
when there are non-zero weak phases. 

Since charm decays proceed in an environment of virulent final state interactions, an 
absence of strong phase shifts in D decays is the least of our concern, since it would happen 
only 'accidentally'. Yet in the presence of hadronic resonances it becomes even a 'mute' 
point, since the resonance provides a mass dependent strong phase that is predictable in 
most cases and thus actually helps to validate a signal. Resonances then create the more 
favourable scenario. 

Ideally we would apply the method proposed by us (see below) to real primary data. 
Unfortunately we do not have access to those. Therefore we start out by using models that 
are consistent with existing data to create a Monte Carlo Dalitz plot; for K^-k^ti^ 
we have been thus 'inspired' by Babar [6] and for — >• tt^tt+tt^ by E791[14]. Then we 
create by hand a single 'seed' for a CP asymmetry and analyze whether our method can 
uncover it; subsequently we vary that single seed. 

We will employ the isobar model [15] for constructing Dalitz plots. The amplitude for 
resonant sub-processes is expressed through Breit-Wigner functions multiplied by angu- 
lar distributions as determined by their spins. The amplitudes of all contributing sub- 
processes are combined coherently with complex coefficients. The latter represent free 
parameters that are fixed from the data using a maximum likelihood fit: the magnitudes 
of the complex coefficients arc related to the fractional contributions of each sub-channel 
and their relative phases reflect the flnal state interactions between the resonances and the 
'bachelor' particles. In Eqs. (6,7) we have exemplified the general procedure by writing 
down the amplitude for / f^. These relative phases are treated as constant, 

since they depend only on the total mass of the system, which in this case amounts to the 
mass of the decaying heavy meson. The non-resonant three-body contribution is usually 
assumed to be fiat over the Dalitz plot or at least described by a smooth distribution. 



2.3 The Novel Proposal 

The challenge we have to deal with in comparing Dalitz plot populations is one of unbiased 
pattern recognition. It is thus analogous to one faced routinely by astronomers: they often 
search for something they do not quite know what it is - at least initially - at a priori 
unknown locations and having to deal with background sources that are all too often not 
really understood. This sounds hke a hopeless proposition, yet astronomers have been 
successful in overcoming these odds. Thus we should be eager to learn from them. 

The Pierre Auger observatory has already adopted the same method for statistical 
weighting in their searches for cosmic ray sources, and we propose to follow suit in defining 
a search strategy for CP asymmetries in Dalitz analyses: rather than study the customary 
fractional asymmetry 

^ iV(0 - iV(0 
N{i) + N{i) 
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in particle vs. anti-particle populations N{i) and N{i) for each bin i, respectively, one 
should analyze the significance 




N{i) - N{i) 



(9) 



N{i)+N{i) 



which amounts to a standard deviation for a Poissonian distribution ^. We will demon- 
strate below through Monte Carlo studies of D and B decays that analysis of the sig- 
nificance a provides a more robust probe of CP symmetry. We will illustrate how the 
observable ^^Scp is highly effective in filtering out genuine asymmetries from statistical 
fluctuations. 

A final technical comment concerning binning size: in the studies presented below we 
have required bins to contain at least twenty events. This number appears 'reasonable', 
but is somewhat ad-hoc. Applying our method to real primary data in the future should 
shed light on the appropriateness of this lower bound. 

2.4 First Summary of the Advantages of Our Proposal 

Analyses of Dalitz plots have so far not 'bagged' any success in establishing CP violation. 
Even so we expect them to become central probes of CP invariance due to the following 



• Local asymmetries are bound to be larger than integrated ones thus facilitating the 
task of controlling systematic uncertainties. 

• The latter - either due to production asymmetries or to detection inefficiencies - 
can be probed and controlled through the analysis of ratios of particle yields. 

• The bin observable ^^Scp defined in Eq.(9) does not suffer from any model depen- 
dance and allows a robust search for asymmetries that are small or in relatively 
small samples. 

• This procedure does not represent a diversion on the (long) path to the ultimate 
goal, namely to arrive at a complete Dalitz plot description and all it can teach us. 
On the contrary - it will accelerate our progress on that journey providing us with 
increasingly powerful pointers for where to focus our attention and constraints for 
the Dalitz parametrizations. 

In the following we will present case studies of B and then D decays to illustrate the 
general method. 

^We will refer to analyzing ^^^Scp instead of A{i) as adopting the 'Miranda' procedure or as 'miran- 
dizing' the CP search. 



features: 
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3 B Decays 



3.1 General Remarks 

Decays B — > /ii/i2^3 with hi = n, K exhibit a pattern in their Dahtz plots that at first 
sight might look surprising: the bands near the edges are crowded while the interior is 
sparsely populated. Yet on second thought this is as expected. For the phase space 
available in B decays is quite large, in particular for non-charm final states. Those will 
typically consist of significantly more than three stable mesons. For three meson final 
states the two primary qq clusters produced in the B decays have to recede from each 
other quickly with untypically low masses; thus they generate the pattern sketched above. 

That final states consisting of just two or three pseudoscalar mesons are a rather 
untypical subset of nonleptonic B decays can be seen also in another way: it has been 
firmly estabhshed that the lifetime of charged B mesons exceeds that of neutral ones [16] : 
t{B~^)/t{B,i) = 1.071 ± 0.009 — in agreement with already the first ftilly inclusive the- 
oretical treatment based on the operator product expansion, which traces this difference 
back mainly to a destructive interference in nonleptonic B^ decays [17]. Yet when one 
sums over the B — > Dn, B — > D*7r and B — > Dp channels one finds that there the B^ 
width exceeds that for B^ by about a factor of two! This is in marked contrast to the 
case of D mesons where the sum of the partial widths for D — > Kit, K*Tr and Kp already 
exhibit the same pattern as the total widths. 

There are many modes that carry considerable promise to reveal CP violation and 
shed light on the underlying dynamics. We will focus on just one B (and later on just one 
D) mode in this note for two reasons: the pedagogical one that we do not want to 'over- 
feed' the reader; and the very practical one that so far little experimental information 
exists about these B decays. In this spirit we will discuss B^ JC^vr+vr". This channel 
is predicted to have a large component from the Penguin operator. Since that operator 
is derived from a loop process — i.e., a pure quantum effect — it represents a wide 
gateway for New Physics. One should also note that the neutral two-body counterpart 
Bd — > K^7r~ has already shown a direct CP asymmetry [18]. 

Observing a CP asymmetry here is unlikely to be the main challenge - that role is 
reserved for the question whether an observed signal is generated by CKM forces alone 
or requires the intervention of New Physics that probably provides merely a non-leading 
contribution. We know of no model-independent way to settle this issue and thus have to 
rely on theoretical treatments that are based on more than just basic features of QCD, 
yet still require model assumptions not (yet) derived from QCD. 

3.2 K^TT+TT- 

Wc will describe this case in considerable detail, since it commands a relatively large 
branching ratio compared to other charmless final states and there is strong evidence for 
a direct CP asymmetry associated with the 5^ — > ii'^p°(770) sub-channel [6, 7]. It also 
provides a clear illustration of the power of our method. 
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The moduli and phases of its amphtudes are 'inspired' by BaBar's results [6]. We 
include five resonant and one non-resonant contribution; the latter is assumed to be fiat 
over the Dalitz plot purely for reasons of convenience and the lack of a specific alternative. 
We analyze two versions each with a single seed of CP violation, namely one with a 
CP asymmetry in the overall phase for the p°(770) and the other one for the /o(980). 
To provide a clear demonstration of our method we start out by assuming the phase 
of the 7^+p°(770) relative to B~ 7^~p°(770) to be large, namely 60", which 

is still allowed by the data [6]. Then we analyze two cases with a significantly smaller 
phase difference, namely 20° and 10°, respectively. In the latter two cases neither a visual 
inspection of the Dalitz plot nor using the fractional difference A{i) suffice to estabhsh 
the resulting CP asymmetry. Yet an analysis of the significance ^^Scp allows even to 
locate the origin of the asymmetry in the Dalitz plot. 

3.2.1 Model 

The specifics of this version are listed in Table 1. For diagnostic clarity we pick two sets 



mode 


a+ 


6+ 


a 


6- 


/s:*(890)7r 


1.0 


0.00 


1.00 


0.00 


K{U30)7r 


2.1 


6 


2.1 


6 


p{770)K 


0.9 


-34 


0.9 


26 




1.0 


132 


1.0 


132 


XcK 


0.3 


-143 


0.3 


-143 


NR 


0.6 


-109 


0.6 


-109 



Table 1: Magnitudes and phases, in degrees, of the amplitudes defining Model " p^" for our 
toy Monte Carlo sample. The difference in the p(770) phase for the S+ and B" channels 
provides the only source for a genuine CP violation. 

of amplitudes for 5+ and B" decays, shown in Table 1, that differ in a single parameter 
only, namely the phase of the p{770)K amplitude, while all moduli of the amplitudes are 
the same. 

With these parameters the signal amplitudes for the B~^ and B~ are integrated over 
the Dalitz plot, yielding a direct CP asymmetry of about 3 x 10^'^. For a sample with 
300K B^ — > K^n^n~ decays this corresponds to 298K B~ —>■ K'n'^n' events. At first 
one might think that the fully integrated rate can show no difference for the S+ and B~ 
channels, when the only seed of CP violation planted into the Monte Carlo model is a 
difference in the overall phase of the p contribution. Yet the small direct CP asymmetry 
is due to the interference of the triangle - pun intended - of Kp, KJq and K*n amplitudes. 

For the B^ and B^ samples we assume a background of about 200K events. The 
resulting Dahtz plots are shown in Fig.l. They do not look quite the same. To make 
their differences more explicit we have plotted the fractional asymmetry A(i) of Eq.8 bin 
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Figure 1: Dalitz plot distribution for S+ — > K'^n'^n (top) and B K tt+tt (bottom) 
in Model 

for bin in Fig. 2. The resulting display is a very noisy one with many bins showing sizable 
differences, both in the p — fo interference region, where our model has to yield a genuine 
asymmetry, and in the central region, where it cannot 

The 'eager' eye might notice that the differences in the former follow a slightly more 
systematic pattern than in the latter, yet it could not be called compelling, in particular 
if we did not know the underlying dynamical structure. 

The effect of the statistical fluctuations can be ilustrated by the following exercise. 
We plot in Fig. 3 the significance distribution for a situation where the 5+ and B~ Dalitz 
plots were generated with exactly the same set of parameters. In this case only statistical 
fluctuations arc observed. The upper display shows most bins to exhibit some differences; 
yet the fact that the ^^S^p distribution is completely consistent with a pure Gaussian 
pattern, as shown in the lower display, reveals them to be consistent with mere statistical 
fluctuations. 

After our method has successfully passed this null test we return to the model deflned 

by Table 1 and Fig 1. To obtain a clearer picture we 'mirandize' our analysis, i.e. turn to 
the significance ^^S^p defined in Eq.9. We plot the resulting values for ^^Scp in the upper 
display of Fig. 4 and its distribution together with a Gaussian fit in the lower display. 
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(B+ - B-)/(B+ + B-) 



Figure 2: Asymmetry in the Dalitz plot bins for Model as defined by Table 1. 

The Dalitz plot of the significance ^^S^p shows a considerably less noisy pattern than 
before with an obvious asymmetry surfacing in the p - fo interference region. The fact 
that a genuine CP asymmetry has surfaced in the Dalitz plot is demonstrated in the 
lower display: there is no acceptable Gaussian fit to the ^^S^p distribution meaning the 
asymmetries are considerable larger than can be generated by statistical fluctuations. 

Our ambition has of course to go further than just knowing that somewhere in the 
Dalitz plot there is a true CP asymmetry - we want to determine in which subdomain(s) 
it resides and whether it is due to an interference between neighbouring resonances or 
due to different widths of two CP conjugate resonances. For that diagnosis we divide the 
Dalitz plot region into subdomains. The choice of these subdomains has to be informed 
by our understanding of the significant subprocesses. In the case of — * K^7[~^n^ we 
divide it into the four regions shown in Fig.5: 1 and 11 containing the p(770) resonance, 
III with the Kit resonances and IV populated mainly by background. In Fig. 6 we have 
plotted the ^^Sc^p distributions separately for these sub-domains. The results are very 
telling: the plot clearly reveals that the asymmetry resides in regions I and 11, while 111 
and IV show no trace of a genuine CP asymmetry - in full agreement with the underlying 
model chosen to generate these Dalitz plots. 

3.2.2 Model "/o" 

In this version we use the same model parameters as above (see Table 1) with two essential 
differences: S+ ^ ^ -34" for the p°(770); S+ = 132" 7^ 5" = 69" for the /o(980), 
i.e., a phase difference of 63 °. Again such a difference is quite compatible with BaBar's 
findings [6]. 
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(B? - B-)/sqrt(B+ + B-) 





Entries 1987 


3_ 


X^/ndf27.74 / 23 
PI -0,2356E-02± 0.2326E-01 
P2 1.017 ± 0.1628E-01 
,P3 1962. ± 44.33 


2_ 






r , , , , / 


/ \ 

< < < 1 < < < 1 < < < 1 < < < 1 < < < 1 < < < 





-5 -4 -3 -2 -1 1 2 3 4 5 

^CP 

Figure 3: Top: Significance ^^S^p plot for two CP conserving 300K signal + 200K back- 
ground samples for CP symmetric decays. Bottom: Caussian fit for the ^^S^p distribu- 
tion; PI, P2 and P3 denote the fit values for the central value, width and normalization 
parameter, respectively. 

The 5+ and B~ Dalitz plots are shown in Fig. 7. One can see that the two plots are 
different. Turning to a plot of the fractional asymmetry A{i) shows there are many bin-by- 
bin asymmetries, yet those exhibit again a rather noise pattern, sec Fig. 8a. Once again 
'mirandizing' the display, i.e., plotting ^^Scp instead of A(i), leads to a more organized 
message, shown in the upper display in Fig. 8b. In particular, when looking at the ^^Scp 
distribution of Fig. 9 we see that over and above the statistical fluctuations there is a 
genuine CP asymmetry. 

As before its location can be narrowed down further by dividing the Dalitz plot in 
the four regions of Fig. 5 and plotting the ^^S^p distributions separately for them, see 
Fig. 9. It clearly identifies regions I and II as the main origin of the asymmetry. That is 
as it has to be, since the interference between the Kp and Kfo amplitudes, which is the 
"engine" of CP violation in our model, takes place mainly there. 



13 



(B?- - B-)/sqrt(B+ + B-) 




Figure 4: Top: Significance °pScp plot for X^tt^tt^ for model "p°". Bottom: 

^^Scp for the bins in Top Figure that pass the statistical cut, fit to a centred Gaussian 
with unit width. PI is the normalization parameter. 

3.2.3 Comparing the "p°" and "/o" Models 

The preceding discussion has shown that the ^pS(;p observable and its distribution pro- 
vides a powerful tool that in a model independent way allows to establish the existence 
of a genuine CP asymmetry over and above statistical fluctuations and even determine 
the subregion(s) of the Dahtz plot, where it originates. For both the two Dalitz models 
employed above it was mainly the p — fo interference domain. 

In addition, a closer analysis allows to distinguish the cases where the asymmetry is 
driven by a difference in the Kp and in the Kfo phase, respectively, for the and B^ 
decays, see Figs. 6 and 9b. The discriminator is provided by the interference with the 
'silent' partner, the K*Tr amplitude. This ability would provide important diagnostics 
about the underlying dynamics: for it would enable us to decide whether the CP odd 
operator generating the asymmetry carries vector or scalar quantum numbers. 
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Figure 5: B"^ K'^tt^tt'^ Dalitz plot for model divided into regions. 




Figure 6: Distribution of figure 4 divided in the regions shown in figure 5. PI is the 
normalization parameter. 
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Figure 7: Dalitz plot distributions for S+ — > K'^n'^n (top) and B K n'^n (bot- 
tom) in Model "/o"- 

3.2.4 The case of a smaller p phase difference of 20" and 10° 

The overall phase differences we have assumed for the two models employed above were 
rather large, although still consistent with present data. Consequently even an unsophisti- 
cated 'look' at the conjugate Dalitz plots suggested the existence of a true CP asymmetry. 
Yet for smaller and presumably more realistic values of these phase differences one needs 
the more refined analysis outhned above. Instead of the 60° phase difference in the p° 
amplitude we had assumed above in our model "p°") we now assume a phase difference 
of just 20° and 10°, respectively, while leaving the other parameters as listed in Table 1. 
Figs. 10 and 11 show the resulting Dalitz plots separately for the and B^ decays; 
they look very much the same now. In Figs. 12 and 13 the lower and upper displays 
show the A(i) and ^^Scp plots, for the two scenarios of a 20° and 10° phase difference, 
respectively. The lower displays of A(i) are very noisy with no clear message. The upper 
display of ^''Scp shows a systematic deviation from zero for the 20° case, while that can 
be hardly said for the 10° case. The existence of a genuine asymmetry is demonstrated 
by the ^^Scp distribution of Fig. 14a. 

Even better, one can localize the region of origin for the CP asymmetry as the one 
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(B+ - B-)/sqrt(B+ + B-) 

Figure 8: Top: Asymmetry in the Dalitz plot bins for Model "/o"- Bottom: Plot of the 
significance ^^Scp for K^n^n^. 

where and /q interference takes place, see Fig. 14b. 

The situation for the 10° case is more delicate. The ^'^Scp plot in Fig. 15a shows there 
is no good Gaussian fit to it: the distribution is a bit wider than a Gaussian expression 
can yield, but still symmetric around its maximum. Yet plotting the ^^Scp distributions 
separately for the four regions as before - Fig. 15b - reveals a clear message: there is a 
true CP asymmetry in regions I and II where p°-/o interference takes place, but none in 
regions III and IV. 

We want to stress that for the last two scenarios - a small phase difference of 20° and 
10°, respectively, the sophistication provided by an analysis in terms of the significance 
^PScp was essential in revealing the underlying dynamics. 

3.3 Future B Studies 

We already mentioned there are several other modes that can be studied with high statis- 
tics by LHCb: 

• — > 7r='=7r+7r~: like B^ K^n'^n' it receives contributions from tree as well as 
Penguin operators, yet with the weight of the former enhanced. It thus represents 
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Figure 9: Top row: Sep for the bins in Fig. 8b that pass the statistical cut, fit to 
a centred Gaussian with unit width for model "/o"- PI is the normalization parameter. 
Bottom two rows: Distribution of top row divided into the regions shown in Fig. 5. PI 
is the normalization parameter. 

a nicely complementary process. 

• The more unconventional channels B"^ rr'^pp, K^pp : the presence of the me- 
son allows us to measure the proton and anti-proton polarization, probing for a 
CP asymmetry, otherwise impossible in two-body decays like — > pp. 

• Bd — Bd oscillations would lead to Dalitz plots for B^ — > Ks7i~^7i~ , where the weight 
of different components would shift with the time of decay thus producing time 
dependent Dalitz plots. 

• The same will happen for Bg — > KsK~7r~^, KsK~^K~ , albeit with a much faster 
oscillation rate. 

We will address these transitions in future work. 

In this note we have shown how mirandizing the analysis of Dalitz plots - i.e., studying 
the 'significance' distributions - can act as a powerful filter against statistical fluctuations. 
Yet real data are also vulnerable to systematic experimental uncertainties. For a full 
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Figure 10: Dalitz plot distributions for K'^n'^Tr (top) and B K n'^n 

(bottom) in a model " p°" with a 20° phase difference. 

validation of our method someone has to apply it to real primary data, to which we have 
at present no access. 



4 D Decays 

The SM generates a relatively dull weak phenomenology for charm transitions: 'slow' 
]jo _ £)Q oscillations and tiny CP asymmetries; this, however, makes it a promising 
landscape to search for New Physics [19, 20, 21]. At the same time wc have to analyze more 
closely how slow is 'slow' quantitatively and how tiny is 'tiny'. One has to concede that SM 
dynamics might saturate the observed size of xd = ^Md/Vu and yo = ^^d/^Di and 
that CKM forces can produce CP asymmetries on the C(10~^) level in singly Cabibbo 
suppressed (SCS) modes. Furthermore, ignoring — oscillations, purely Cabibbo 
allowed (CA) and doubly suppressed (DCS) channels (i.e. those without a Ks or Kl) 
cannot exhibit direct CP violation. Any such effect in DCS modes and one on the about 
0.01 or larger level in SCS decays will thus establish the intervention of New Physics. 
Basing such claims on somewhat smaller effects will require theoretical progress that 
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Figure 11: Dalitz plot distributions for K'^n'^Tr (top) and B K n'^n 

(bottom) in a model with a 10° phase difference. 

appears quite feasible. 

The phase space available in D decays is significantly smaller than in B decays. Two- 
and quasi-two-body channels make up more than half of the full nonleptonic width. Fur- 
thermore the Dalitz plots are populated more thorougly than for B decays. Determining 
the impact of individual contributions therefore amounts to a more delicate task. 

4.1 D+ TT+TT+TT- 

In the SM there are already two different amphtudes contributing to these SCS transi- 
tions, and they carry a relative weak phase, albeit a tiny one ~ O(A^) ~ 10~^. Finding 
CP asymmetries significantly larger than 10~^ would provide strong prima facie evidence 
for the presence of New Physics. Searching down to asymmetries as small as 10^"^ re- 
quires huge statistics as well as excellent control over systematics; the method proposed 
by us should be a powerful tool in taking up such a challenge. Probing the whole Dalitz 
plot with its various structures should allow us to make the case for New Physics even 
compelling making use also of the anticipated theoretical refinements sketched below. 
It should be noted that New Physics scenarios like the Littlest Higgs Model with 
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(B+ - B-)/sqrt(B+ + B-) 

Figure 12: Plot of A(i) (top) and °pScp (bottom) for B+ K+n+n' and B' 
K~7r'^7r~ in a model with a 20" phase difference. 

T parity could have an observable impact here through new heavy states appearing as 
virtual particles in Penguin diagrams [22]. 

We adopt a decay model containing four components, namely 

• L>± ^ a^n^. 
. L>± ^ /o7r± 

• a uniform non-resonant D"^ — > tt+tt"?!^. 

We have used the values obtained by Fermilab experiment E791 [14] for the magnitudes 
and phases of the amplitude coefficients. 

A difference in the a phase - a very conceivable scenario - will affect many parts of 
the Dalitz plot and induce CP asymmetries, since the a possesses a very sizable width 
relative to the phase space available in D decays. The resulting complexities are very 
intriguing and will be analyzed in a separate paper. 
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(B+ - B-)/sqrt(B+ + B-) 

Figure 13: Plot of A(i) (top) and °pScp (bottom) for B+ K+n+n' and B' 
K~7r'^7r~ in a model with a 10" phase difference. 

The case of the /o amplitude having a different phase in and decays is very 
interesting for another reason: as long as it has any uu or dd component it will contribute 
in this final state. We have found that this relatively small contribution can still produce 
a clear signature in CP asymmetries mainly due to the narrow width of the /q. Details 
of the required analysis will also be given in the future paper. 

In this pilot study we will focus on one scenarios, which leads to clean signatures, 
namely those with a 1% (equivalent to 3.6°) phase difference in p° amplitude. We have 
selected much smaller phase differences here than for our discussion of B decays above 
for two reasons: (i) They represent much more realistic New Physics scenarios, (ii) Due 
to the considerably smaller phase space and thus shrunk Dalitz plot areas one can expect 
such effects to be still observable. 

In Fig. 16 we display the Dalitz plot for this model and the ^^S^p distribution for 
the whole plot as well as for the two regions I and II. The overall ^^Scp distribution 
unequivocally reveals the existence of a CP asymmetry, since it does not at all follow 
a Gaussian fit. The distributions for the two regions I and II exhibit a very telling 
pattern, namely a rather asymmetric distribution of the bin-wise CP asymmetries. That 
is how it has to be for the line dividing I and II chosen to go through the gap between 



22 




Figure 14: Top row: Distribution of ^^Scp that pass the statistical cut, fit to a centered 
Gaussian with unit width; PI is the normahzation parameter. Bottom two rows: Dis- 
tribution of °^Scp divided into the regions shown in Fig. 5 in a model with a 20" 
phase difference. PI is the normalization parameter. 

the two p 'lobes', as can be seen by straightforward arithmetic. Applying Eqs. (6,7) 
to — > 7r^7r''"7r~, with an analogous, though smaller p - fo interference used in those 
equations leads to following difference in the D'^ — > tt+tt+tt" - D~ — > 7r~7r''"7r~ amplitude 
squared: 



AM = |M+|^-|M-P = [(a^f -(a^_f]|F^^^|^cos2^ + [(a(f -(a^^llFf^l^ 

+2cOS^|Fj^W|2|^BW|2^ 

{[{^1 - s){m} - s) - mprpmfrf]y+ai cos(5^ - Si) - a'_aL cos{S'_ - Sl)] 
-[mprp{m} -s)- mfFfiml - s))K«i sin(5^ - 5^) - a^_aL sm{S'_ - 6L)]} 

(10) 

The term quadratic in cos 6 is responsible for the parabolic shape of the spin one reso- 
nance seen in Fig.l. Yet the interference generates a term linear in cos^. Therefore the 
interference is destructive in region I of Fig. 16- thus implying fewer events for than 
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Figure 15: Top row: Distribution of ^^Scp that pass the statistical cut, fit to a centred 

Gaussian with unit width; PI is the normahzation parameter. Bottom two rows: Dis- 
tribution of °^Scp divided into the regions shown in Fig. 5 in a model with a 10" 
phase difference. PI is the normalization parameter. 

D~ - and the opposite in region II. 

This example illustrates the power of the mirandizing procedure to unequivocally 
uncover even a small asymmetry and track its local origin in the Dalitz plot. 

4.2 Future D studies 

As before with B decays many promising channels await careful study: 

• The more complex scenarios in — > n^n'^n'', where the seeds for CP violation 
reside in the /o and a amplitudes deserve detailed analysis. 

• The doubly Cabibbo suppressed modes — > K^7r'^7r~, K^K^K~ could reveal a 
new source of direct CP violation [23] . 

• With the observation of — oscillations one expects time dependent Dalitz 
plots to emerge in D° — > Kst^'^t^" ■ This time evolution will allow to differentiate 
between direct and indirect CP asymmetries. 
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Figure 16: Dalitz plot for D"^ -K^-K^-n~ in a model with a 1% (3.6°) phase difference 
in the p° amplitude with sub-domains I and II; distributions of the significance ^^Scp for 
the whole plot and the two sub-domains I and II; PI is the normalization parameter. 



5 On Refining the Theoretical Tools 

Even lattice QCD does not allow to treat final state interactions as a matter of principle, 
except for kaon decays, where elastic unitarity can be assumed. Elastic unitarity makes 
little sense for B decays; it might be an approximation of some value in D decays, but we 
have no reliable even semi-quantitative estimate for how good an approximation it might 
be. 

We should be able to clarify the picture at least somewhat by adopting 'theoretical 
engineering' [24]: One considers — > PP {P = pseudoscalar meson) on all Cabibbo 
levels for D'^, and mesons. Relying on a modicum of theoretical judgement 
one selects diagrams deemed relevant for these processes and expresses their amplitudes 
in terms of the known CKM factors and radiative QCD corrections and the a priori 
unknown moduli and strong phases of their matrix elements. Fitting these expressions 
to a comprehensive body of well measured branching ratios one fixes these moduli and 
strong phases. The resulting owe?x;onstraints provide a check on the reliability of such a 
fit. The analogous procedure is then applied to D(^s) PV (V = pseudoscalar meson). 
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While such an analysis cannot replace a full Dalitz plot description, it can provide valuable 

constraints on the latter. 

Alternatively we should be able to develop some framework where we can have a semi- 
quantitative treatment of the interference of a narrow resonance with a broad non-resonant 
contribution that to first approximation can be considered even as flat. 

6 Summary and Outlook 

So far Dalitz plot studies have not established any CP violation with at least five sigma 
significance — yet we are confident this period will soon come to an end. We actually 
expect such studies to become a central tool for obtaining a more detailed picture of 
and perspective on limitations of CP invariance. An acceptable description of the Dalitz 
plot usually has to satisfy a sizable number of owerconstraints, which provides a powerful 
validation tool to control systematics. Furthermore — and maybe even more importantly 
— it provides us with information about the Lorentz structure of the underlying transition 
operator that cannot be inferred from partial rate asymmetries in two-body final states. 

A full fledged Dahtz plot description thus represents the 'holy grail' in our quest for 
mapping out CP violation in B, D and maybe even top quark decays. The journey there 
will however require a substantial amount of time, as it is with all 'holy grails'. It also re- 
mains to be seen to which degree there will arise uncertainties due to an irreducible model 
dependance. The method we have proposed in this note for searching for CP asymmetries 
in the populations of Dalitz plots is not meant to replace Dalitz plot parametrizations: 

• The proposed method will allow to establish the existence of CP asymmetries with 
more limited statistics and identify their topography in the Dalitz plot in a robust 
and model independent way. 

• Furthermore isospin sum rules [25] can already be applied to its findings. 

• It will speed up the construction of the full Dalitz plot description and provide 
powerful validation for it. 

In this paper we have described a model independent method for establishing the 
existence of a CP asymmetry in a Dalitz plot and inferring its location. To fully gauge its 
power it is important to apply it to high statistics primary data with their experimental 
systematic uncertainties. If it passes that test, then one can study how to extract maximal 
information about its parameters, in particular the weak phase producing it. The relevant 
expression is given in Eq.(lO). Various methods can be employed to achieve such a goal; 
finding the optimal one requires future detailed analysis. ^ Various features can be 
employed to discriminate between SM and potential New Physics effects. In B ^ Ktiti 
the SM can generate a significant weak phase only through its {V — A)x {V — A) currents, 
since the 6 — > s Penguin operator does not carry a weak phase. New Physics thus could 
make its ]:)resence felt through ]:)roducing a weak phase for a scalar state like the /q. 

^We thank J. Appel for forcing upon us an illuminating discussion of this point. 



26 



Clearly a large amount of also theoretical work is required. While we should not count 
on theorists achieving miracles, we can expect a positive learning curve for them. 
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